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Abstract. This chapter reviews the following items: 

1. Energy deposition and track structure of fast positrons: ionization slowing down, 
^ ' number of ion-electron pairs, typical sizes, thermalization, electrostatic interaction 
q . between e + and its blob, effect of local heating; 

2. Positronium formation in condensed media: the Ore model, quasifree Ps state, 
intratrack mechanism of Ps formation; 

^- . 3. Fast intratrack diffusion-controlled reactions: Ps oxidation and ortho-para 

conversion by radiolytic products, reaction rate constants, interpretation of the PAL 
-Y~. ■ spectra in water at different temperatures; 

l/~) , 4. Ps bubble models. "Non-point" positronium: wave function, energy 

~f~, ' contributions, relationship between the pick-off annihilation rate and the bubble radius. 
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Positrons (e + ) as well as positronium atoms (Ps) are recognized as nanoscale probes 
of the local structure in a condensed phase (liquid or solid) and of the early radiolytic 
physicochemical processes occurring therein. The parameters of positron annihilation 
spectra determined experimentally (e.g., positron and Ps lifetimes, angular and energetic 
widths of the spectra, Ps formation probability) are highly sensitive to the chemical 
composition, the local molecular environment of Ps (free volume size), and the presence 
of structural defects. They are also sensitive to variation of temperature, pressure, 
external electric and magnetic fields, phase transitions. 

The informative potentiality of positron spectroscopy strongly depends on the 
reliability of any theory describing the behavior of positrons in matter, since it 
should help to decipher the information coded in the annihilation spectra. So, 
realistic (quantum mechanical) models are needed for e + track structure, e + energy 
losses, ionization slowing down and thermalization, intratrack reactions (ion-electron 
recombination, solvation, interaction with scavengers), Ps formation process, Ps 
interaction with chemically active radiolytic species and e + /Ps trapping by structural 
defects. 

Usually, treatment of the measured annihilation spectra is reduced to their 
resolution into a set of simple trial functions: sums of decaying time exponentials in 
the case of PALS (Positron Annihilation Lifetime Spectroscopy), of Gaussians in the 
case of ACAR (Angular Correlation of Annihilation Radiation) and DBAR (Doppler 
Broadening of Annihilation Radiation). The outcome of such conventional analyses of 
positron annihilation data are the intensities of these components and the corresponding 
lifetimes/widths. 

However, realistic theoretical models suggest more complex kinetics for describing 
the intratrack reactions between the positron and Ps atom and radiolytic products 
and scavengers. For example, because of the inhomogeneous spatial distribution of 
the track species, their outdiffusion becomes an important factor: diffusion kinetics 
cannot be expressed in terms of mere exponentials or Gaussians. Obviously, elaborated 
theoretical models should be used in the fitting procedure of the annihilation spectra. 
In this case, the adjustable/fitting parameters (reaction rate constants or reaction radii, 
diffusion coefficients, initial size of the terminal part of the e + track, contact density in 
the Ps atom) would present a clear physical meaning instead of the above mentioned 
"intensities" of some trial functions. 

The present chapter is an attempt to a systematic presentation of the theoretical 
model, which describes e + fate since its injection into a liquid until its annihilation. It 
should be a basis for a positron spectroscopy of liquid media. 
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Figure 1. Scheme of the end part of the e + * track and Ps formation. 



1. Energy deposition and track structure of the fast positron 



1.1. Ionization slowing down 

Positrons, produced in nuclear /3 + -decay, have initial energies of about several hundreds 
of keV. Once injected into a medium they lose energy via molecular ionization. Within 
10 ps the positron energy drops down to the ionization threshold. Further approach 
to thermal equilibrium proceeds primarily via excitations of intra- and intermolecular 
vibrations. This usually takes several tens of picoseconds [TJ [2] . 

Roughly half of the positron kinetic energy is lost in rare head-on collisions, resulting 
in the knocking out of ^-electrons with kinetic energies of about several keV, the tracks 
of these electrons forming branches around the positron trail (Fig. [I]). The other half of 
the energy is spent in numerous glancing collisions with molecules. The average energy 
loss in such a collision is 30-50 eV (at maximum, 100 eV). A secondary electron knocked 
out in a glancing collision produces, in turn, a few ion-electron pairs inside a spherical 
microvolume, called a "spur" in radiation chemistry. Its radius, a sp , is determined by 
the thermalization length of the knocked-out electrons in the presence of the Coulombic 
attraction of the parent ions. Note: why restrict to the electrons knocked-out in a 
glancing collision? 

While the positron energy W is greater than W cy i ~ 3 keV, the mean distance k 
between adjacent ionizations produced by the positron is greater than the spur size 2a sp . 
It ensues that at high positron energies, the spurs are well separated from each other. 
The trajectory of the fast positron is a quasi-straight line because U is less than the 
positron transport path l tr . The latter is the mean distance traveled by the positron 
before it changes its initial direction of motion by 90°. 

When li < 2a sp < l tr or Wu < W < W cy i the spurs overlap, forming something 
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like a cylindrical ionization column. When the e + energy becomes less than the blob 
formation energy, Wu (~ 500 eV, see below), the positron is about to create a terminal 
blob. The diffusion motion of e + in the blob becomes more pronounced: the direction 
of its momentum changes frequently due to elastic scattering and the ionization of 
surrounding molecules. Roughly speaking, all intrablob ionizations are confined within 
a sphere of radius au- The terminal positron blob contains a few tens of ion-electron 
pairs (no ~ Wu/Wi ep ~ 30) because the average energy Wi ep required to produce one 
ion-electron pair is 16-22 eV [3]. Finally, the positron becomes subionizing and its 
energy loss rate drops by almost 2 orders of magnitude |lj. 

Typical dependencies of k(W) and k r (W) versus the energy of the positron are 
shown in Fig. [2l The calculation of k(W) = Wi ep /LET(W) is based on LET (Linear 
Energy Transfer) data of e 1 * 1 . The estimation of the transport path has been done in 
the framework of the Born approximation (the wavelength of the positron with energy 
> 100 eV is small in comparison with the size of molecules), where the Born amplitude 
was calculated by simulating a molecule of the liquid as an iso-electronic atom. 

At low energies, positron scattering becomes more and more efficient and its motion 
must be regarded as diffusion- like with the energy- dependent mean free path, l tr (W), 
between successive "collisions" , which completely randomize the direction of the velocity 
of the particle. During ionization slowing down, when e + * loses energy from Wi down 
to Wf the mean square displacement of the positron is 
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where D P (W) = l tr (W)v p /6 is the energy- dependent diffusion coefficient and dt = ^ 
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v p \—dW/dx\i . 

following equations: 
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Now we can define Wu and au- These quantities are obtained from the 
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Here Ry=13.6 eV stands for a typical ionization potential. Eq. ([2]) indicates that the 
terminal positron blob is a spherical nanovolume, which confines the end part of its 
trajectory. This is where ionization slowing down is most efficient (the thermalization 
stage of the subionizing positron is not included here). The mathematical formulation of 
this statement is twofold. Just after the first "step" of the blob formation, which is equal 
to ltr(Wu) (the thick red arrow in Fig. [TJ, the positron reaches the center of the blob. 
After that, the end part of the ionization slowing-down trajectory is embraced by the 
blob; i.e., the mean-square slowing-down displacement of the positron, R 2 on (Wu, Ry)— a| 
is equal to a 2 t the "radius" of the blob squared. 

The solution of these equations in the case of liquid water is shown in Fig. [2j It 
gives Wu ~ 500 eV and a bL w 40 A. One may assume that the values of a bL and Wu do 
not differ significantly from one liquid to another because the ionization slowing-down 
parameters depend mostly on the ionization potential and the average electron density, 
parameters that are more or less the same in all molecular media. At first sight it 
seems strange that the obtained blob size au ~ 40 A is smaller than the dimensions of 
ion-electron pairs (60-200 A), well known from radiation chemistry [5J. This is because 
the dimension of the ion-electron pair is determined by the thermalization length of the 
subionizing electron (in the field of the parent cation), while the positron blob radius is 
determined by the ionization slowing-down of the energetic positron when it loses the 
last ~ 500 eV of its energy. 

1.2. Thermalization stage. Interaction between the positron and its blob 

At the end of the slowing-down by ionization and electronic excitation, the spatial 
distribution of e + * coincides with the distribution of the blob species (i.e., ~ 
exp(— r 2 /a^)). Such a subionizing positron having some eV of excess kinetic energy 
may easily escape from its blob because the interaction between the blob and the e + * is 
rather weak (see next section). It is expected that by the end of thermalization, the e + 
distribution becomes broader with the dispersion: 

a 2 p ^a 2 u + (R 2 vib (Wo,T)) Wo . (3) 

Here Rl ib (T, W ) is determined by Eq. ([1]), where \dW/dx\ ion should be replaced by 
\dW/dx\ V ib, i.e., the stopping power of subionizing e + * related to the excitation of 
vibrations (T is the temperature in energy units). The estimation of a p requires 
quantitative data on \dW/dx\ vi h, scattering properties of subionizing e + *, and the 
spectrum of its initial energies Wo after the last ionization event. In Eq. ([3]) (. . .)w 
denotes the average over Wq. In contrast to the parameters related to ionization slowing 
down, a p strongly depends on the properties of each particular liquid and may reach 
hundreds of A especially in non-polar media. 

Between the positively charged ions and the knocked-out intrablob electrons there 
exists a strong Coulombic attraction: out-diffusion of the electrons (even during their 
thermalization) is almost completely suppressed and the distribution of the ions is close 
enough to that of electrons. This case is known as ambipolar diffusion when ions and 
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electrons expand with the same diffusion coefficient equal to the duplicated diffusion 
coefficient of the ions. Thus, the distribution of the blob electrons remains slightly 
broader than that of ions. So, the calculated electrostatic potential therein is everywhere 
repulsive towards the positron and its typical value is about several T. 

However, there is an opposite effect: while residing inside the blob, the thermalized 
e + rearranges the intrablob electrons, so that the total energy of the system decreases 
because of the Debye screening. The corresponding energy drop may be estimated 

2 

using the Debye-Huckel theory. It is ~ — 573- . However, the Debye radius 

e(r D +a M /n < ) 

rr, 7n (Aiir c Ci ep )~ l l 2 ~ 4 A is quite small in comparison with au/n the average 
distance between intrablob electrons (here r c = e 2 /e OQ T is the Onsager radius and 
Ciep ~ n o /^ na bi i s the concentration of ion-electron pairs within the blob). Thus, this 
screening energy ~ Tn r c /au becomes a dominant contribution, which is some tenths 
of eV. It may result into e + trapping within the blob. 

The above estimation shows that both the thermalized and epithermal positrons 
could be trapped within the blob. This effect is important to calculate the Ps formation 
probability in an external electric field [HI [7] . 

1.3. Effects of local heating and premelting in the terminal part of the e + track 

We have estimated above that the ionization slowing down of energetic positrons (e + *) 
and subsequent ion-electron recombination release an energy up to 1 keV in the terminal 
e + blob, this energy being finally converted into heat. Therefore the temperature in the 
e + blob should be higher than the bulk temperature. This phenomenon may be called 
as the local heating effect. This transient temperature regime may strongly affect, for 
example, the Ps bubble growth, by changing viscosity of the medium[8]. This effect may 
also influence the mobility of intratrack species and their reaction rates constants. 

For quantitative estimations we simulate this process with the help of the 
macroscopic heat transfer equation: 

dT(r t) 
c P p K m ' = div(AVT) + q + (r, t), T{r, t = 0) = T bulk . (4) 

Here T(r,t) is the local temperature, T^uik is the bulk temperature of the medium, 
c p is its specific heat capacity, p is the density, A is the thermal conductivity. The 
second term in the RHS quantifies the energy released by the positron when creating 
its blob: q + (r,t) ~ WbiG(r,a)f(t,r), where Wu ~ 1 keV is the blob formation 

-7' 2 /a 2 

energy, G(r, a) « e 3/2 3 " describes the spatial distribution of the released energy and 

"" a bl 

f(t, t) ~ exp ( — - ~ r 2 j /(\^2ttt) is its temporal distribution, where r w 0.3 ps is the 
typical time of ion-electron recombination. 

In general, the system can be simultaneously solid and liquid. Then, the following 
method for obtaining a numerical solution of Eq. (J3J) can be used [TU [15] . To describe 
the deposition of the latent heat of melting, q m , we added to c p {T) a term, which is 
non-zero only in a narrow temperature interval T m — AT < T < T m , where T m is the 
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Figure 3. Calculated temperature profiles in ice T^uik = —10 °C (on top) and in 
water Tb u ik (bottom) is slightly above °C. au = 40 A, Wu = 1 keV. 



melting point temperature and AT is the width of the phase transition (arbitrarily fixed 
to 0.25 K). Moreover, an integration of this term over temperature must give q m . This 
additional overshot to c p (T) is simulated by a Gaussian function as follows: 

(T - T m f 



c p (T) = c p (T) +g m exp 



V2nAT 2 



(5) 



2AT 2 
Close to the phase transition region, the T- dependencies of thermal conductivity and 



Table 1. Thermodynamical properties for water and some alcohols in solid and liquid 
phases near melting temperature [9] [TO] [IT] [T2j [13] 



substance 


T K 

- 1 mi - lv 


X s / \ L ; W/m/K [IU 12] 


q m ; J/g [13J 


p S J p L. g / cm 3 ^Qj 


methanol 


175.6 


0.32/0.21 


99 


0.98/0.79 


ethanol 


159 


0.27/0.17 


109 


1.06/0.79 


butanol 


184.5 


0.47/0.16 


125 


1.05/0.81 


water 


273 


2.38/0.56 


334 


0.92/1.00 



Positronium in a liquid phase 



0.8 

" 0.6 

c 

J* 0.4 

0.2 





M^ m 






















o 


butanol 


J* 


c d jto: 




V 


□ 


ethanol 


&rn to< 


ffiK&rO 






V 


methanol 


Wo>/ 










water 







-10 



-6 -4 

T-T ,K 



70 
65 
60 
55 
50 
45 



40 L 
-10 








O 


butanol 


. 




D 


ethanol 






V 


methanol 
water 


, 


, 






6 


-4 


-2 


T-T 


, K 







Figure 4. Lifetime of the molten region t m (left) and its radius R m (right) for 
methanol, ethanol, butanol and water are displayed as a function of X), — T m . 
density were approximated by smooth functions: 

X s X L 
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exp(-(T m - T)/AT) + 1 + exp((T m - T)/AT) + 1 
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P 



exp(-(T m - T)/AT) + 1 exp((T m - T)/AT) + 1" 
This approach works well especially when it is hard to trace an interphase boundary 

EH- 

Simulations were done for water, methanol, ethanol and butanol. Some 
thermodynamic properties of these substances are given in Table [U 

Temperature profiles T(r,t) were calculated numerically for < r < 200 A and 
< t < 1 ns. Eq. (j4j) was solved as a 1-D problem with boundary conditions 
T(r = 200 A) = Tbuik with the help of PDEPE subroutine from Matlab. Some 
temperature profiles in ice and in water close to the melting point are shown in Fig. |3j 

At the transition, the phase of the medium is deduced from the temperature of each 
spatial point, which can be higher than T m + AT (liquid phase) or lower than T m — AT 
(solid phase). Fig. @] displays the maximum radius of the molten region, R ma x, vs. Tb u ik 
and the lifetime t max of the molten region (at t > t max temperature of an any point of 
the medium is below T m + AT). 



2. Positronium formation in condensed media 

2.1. The Ore model 

The Ore model was proposed for interpretation of the Ps formation in gases [16J. It 
implies that the "hot" positron, e + *, having excess kinetic energy, pulls out an electron 
from a molecule M, thereby forming a Ps atom and leaving behind a positively charged 
radical-cation M + ': 

e + * + M^Ps + M + \ (7) 

This process is most effective when the e + * energy W lies within the "Ore gap" : 

I G - Ry/2 < W < W ex (otI g ). (8) 
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Figure 5. Energy diagram of the Ps states. Energy of the Ps atom in vacuum is —6.8 
eV. Vq Ps is the Ps work function (the energy needed for Ps to enter the liquid without 
any molecular rearrangement) . qf-Ps senses the repulsive potential from the cores of the 
molecules (for convenience, this potential is shifted down by 6.8 eV). Being averaged, 
the latter may be represented as a step of height V^ s = U . E^ sw Ry/2e 2 is a rough 
estimation of the qf-Ps energy in a dielectric continuum. Within some picoseconds, 
qf-Ps transforms into the bubble state with the free-volume radius, R, the center-of- 
mass wave function of the Ps, ip(r), and the ground-state energy Ep s — Ry/2 (for large 
bubbles). Rjj is the radius of the potential well, where Ps becomes localized. 



Here Iq is the first ionization potential of the molecule, W ex is its electronic excitation 
threshold and Ry/2 = 6.8 eV is the Ps binding energy in a gas phase. If the positron 
energy is lower than Iq — Ry/2, e + cannot pick up an electron from a molecule. When 
W > W ex , electronic excitations and ionizations dominate and Ps formation becomes 
less effective. 



2.2. Quasifree Ps state 

In the condensed phase, because of the presence of molecules and lack of free space, the 
final Ps state differs from that in vacuum [17J. We call this state quasifree positronium 
(qf-Ps). If we adopt a binding energy Ej, of qf-Ps in a dielectric continuum of roughly 
Ry/2e 2 (instead of just Ry/2 in vacuum), where e = n 2 w 2 high-frequency dielectric 
permittivity is the square of the refractive index, it is seen that the Ore gap in a medium 
gets squeezed 

I L - Ry/2e 2 < W < W e L x (9) 

and may even completely disappear because of a significant decrease in the binding 
energy. Here, W^ x is the lower threshold of electronic excitations of a molecule in a 
liquid, II = Iq — V$ — \U P \ is the liquid phase ionizing potential of a molecule (Vq is 



Positronium in a liquid phase 10 

the electron work function, or the energy of the ground state of electron in a medium 
and U p is the energy of polarization interaction of the positively charged ion with the 
medium) . 

In [18] it was suggested that in some polymers, because of the presence of large 
free-volume space, the final Ps state may not be quasifree: Ps would form immediately 
in one of the preexisting voids. If it did, the total energy gain would be larger and 
the Ore gap could exist. However, because the final Ps state has obviously zero 
translational momentum in contrast with the projectile "hot" positron, this reaction 
must be suppressed owing to the conservation law of momentum. e + * cannot promptly 
transfer its momentum to surrounding molecules because of the large difference in 
masses. Momentum relaxation as well as energy absorption during formation of the 
localized Ps in a cavity cannot be immediate and requires some time. Probably, we 
should expect the formation of qf-Ps in the bulk at first. Then this particle can lose 
energy and momentum before eventually ending up in a suitable cavity after a while. 

2. 3. Recombination mechanism of Ps formation 

The above discussion points out the unique possibility for Ps formation in molecular 
media. This mechanism postulates that Ps is formed through the combination of the 
thermalized particles - quasifree positron and one of the intratrack electrons: 

e qf + e biob -^ (e + ---e") ->• qf-Ps -» Ps in the bubble. (10) 

This reaction proceeds in the terminal part of the e + track (in the e + blob). If the 
positron is thermalized outside the blob, the only way for it to form Ps is to diffuse 
back and pick up one of the intrablob electrons. Otherwise it annihilates as a "free" 
positron. When e + picks up an electron, the "initial" separation between them must be 
comparable to the average distance (47ra^/3n ) 1 / 3 ~ 20 A between intrablob species. 
The total energy of this pair is approximately the sum V^ + + Vq of the e + and e - work 
functions. In comparison with the work functions, the binding energy of such a e + ■ • • e~ 
pair is small (about 0.1 eV). However, the translational kinetic energies of the particles 
must be less than the binding energy, otherwise the pair will break up. So, just before 
Ps formation, the positron and track electrons must be almost thermalized. 

However, such a weekly bound pair is not at the bottom of its energy spectrum, Fig. 
EJ The two particles approach each other (in average) and continue to release energy via 
excitation of molecular vibrations. Finally, the pair reaches an equilibrium state, which 
we term the quasifree Ps. Roughly, the binding energy of qf-Ps is Eb ~ Ry/2e 2 , about 
1 eV. In liquids, further energy gain of the pair is related with the rearrangement of 
molecules and appearance of some additional free space around Ps. e + and e~ get closer, 
repelling molecules from their location and forming a Ps bubble state. A substantial 
decrease in their Coulombic energy is the driving force of this process. 

Over the last 30 years the (re) combination mechanism has become extremely 
widespread [TH [19]. It has been used to interpret numerous data on Ps chemistry, 
and explain variations of the Ps yields (from to 0.7) in very different condensed 
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media, where parameters of the Ore gap are practically the same. It provides a 
natural explanation to the changes in the Ps formation probability at phase transitions. 
Experimentally, the observed monotonic inhibition of Ps yields (practically down to 
zero) in solutions of electron acceptors contradicts the Ore model, but inserts well in the 
recombination mechanism. It explains the anti- inhibition effect, including experiments 
on Ps formation in moderate electric fields in pure liquids and mixtures. 

There are two models which utilize this mechanism, the spur model [20] and the 
blob model (diffusion- recombination model) [2T1 122] . In spite of the fact that both 
models answer the question about the Ps precursor in the same way, they differ as to 
what constitutes the terminal part of the e + track and how to calculate the probability 
of the Ps formation [T7] . 



3. Intratrack reactions in the positron blob 

The interaction of the Ps atom with primary radiolytic products, formed in the terminal 
part of the positron track (in the e + blob) due to e + ionization slowing down, is feature 
inherent to the positron spectroscopy of any molecular medium and must be taken 
into account in interpreting the annihilation spectra. Many intratrack products are 
strong oxidizers and/or radicals (like H30 + radical-cations and OH-radicals in water) 
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and their initial concentration in the e + blob is not small (up to 0.05 M). The average 
distance between intrablob species is comparable with the diffusion displacement of the 
Ps bubble before e + annihilation, so the contribution of diffusion-controlled oxidation 
and ortho-para conversion reactions is quite possible. 

Neglecting these reactions leads to obvious contradictions. Fig. [7] shows the In- 
dependence of the lifetime (r 3 ) of the long-lived component of PAL spectra in pure 
water (roughly, r 3 is the pick-off lifetime of ortho-positronium) . Experimental data 
from different authors reasonably agree, but all are in strong contradiction with the 
theoretical expectations based on any Ps bubble model [23]: with increasing T, the 
surface tension coefficient decreases, so the size of the Ps bubble should increase, which 
should lead to an increase in r 3 with T, in sharp contradiction to the PAL data. 

One most informative way to investigate Ps intratrack reactions is to study the 
temperature variation of the PAL spectra. We have chosen liquid water to initiate this 
investigation, because radiolytic processes in water have been studied for decades and 
are well known. In addition, the Ps bubble growth proceeds very fast in water (< 10 ps) 
and therefore may be considered as an instantaneous process. An opposite situation is 
true in glycerol at T < 50 °C, where Ps interaction with intratrack radiolytic products 
and formation of the Ps bubble state may occur simultaneously jS]. Because all these 
processes strongly depend on temperature, it is worth to consider two extreme cases, 
the cases of "high" and "low" temperatures. 

3.1. High temperature region 

This region may be defined through the following conditions: 

1) the Ps bubble formation time, tb u bbie(T) fa Rrj/a [25] is short as compared to the 
time resolution of a PAL spectrometer, ~ 0.1-0.2 ns. Here, R ~ 3-6 A is the radius of 
the Ps bubble, rj and a are viscosity and surface tension coefficients of a liquid medium. 
Experimental observation of the bubble growth is impossible in this case because its 
equilibrium value is reached too fast (less than 0.1 ns). 

2) the diffusion length a/6(D Ps + Di)r ?J of Ps and radiolysis products (subscript i) 
in the e + blob during the ortho-Ps lifetime, r 3 (about few nanoseconds), must exceed 
the average distance between intrablob particles f m 10-20 A. This can be expressed 
through the condition t di ff < r 3 , where £<&// = r 2 /6(Di + -Dp s ); on the basis of the 
Stokes-Einstein relationship, Di fa Dp s fa kBT/6iTr](T)Ri, Ri is a few A. In glycerol 
both conditions are satisfied at the same temperature region T > 90 °C (Fig. [8]), but in 
liquid water they are fulfilled at any T from to 100 °C . 

In the high-T region the influence of intrablob reactions is important. The radiolytic 
processes initiated through ionizations induced by the slowing down of fast e + * may be 
represented through the following basic reactions: 

ionization: 



A-* 



+ RH^e- + RH + + e + , (11) 
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Figure 8. Temperature dependences 
of the Ps bubble formation time tbubbie 
and typical time tdiff of the diffusion- 
controlled intrablob reactions in glycerol. 



ion-molecule reaction: 

RH+ + RH ->■ RH+ + R, 
ion-electron recombination: 

e" + RH + ->■ -RH* ripl , RH* ingl ; 

e~ + RH+ ->■ RH + H ->• R + H 2 
electronic deexcitation: 

+RH 



RH 



tripl 



2R + H 2 



R H* ingl ->■ RH. 



(12) 

(13) 
(14) 

(15) 



From reactions ( 1Tl~[fT4"|) one can infer that the total number of R in the e + blob varies 
weakly in time and is approximately equal to the initial number, no, of ion-electron pairs 
in the blob. However, their spatial distribution broadens in time due to out- diffusion. 
The main effect of ions and radicals on Ps is oxidation and spin conversion reactions. 
If we denote for simplicity all these species (RH + , RH^ , R, H) by the same symbol, R, 
and term them "radicals" , their reactions with Ps may be written as: 
Ps oxidation: 



Ps + R^e + + R", 

ortho-para conversion: 

1 3 

para- or ortho-Ps + R — > -para-Ps + -ortho-Ps + R. 

4 4 



(16) 



(17) 



In the high temperature region these Ps reactions are mostly diffusion-controlled. 
With the help of the nonhomogeneous chemical kinetics approach (the blob model, 
[2"4"| |2"6] ) reactions ( II llTlT|) may be described in terms of the following equations on the 
concentrations of different intratrack species: 



dcj(r,t) 

dt 



U^LXCi riijC^Cj AiCii 



Ci(r,0) 



exp(-r 2 /qf) 

1 7r 3 / 2 af ' { } 



where Ci(r,t) is the concentration of the intrablob species of the z-th type (including 
all the positron states: e + , para-Ps, ortho-Ps), k^ is the rate constant of the reaction 
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between the i and j reactants and Aj is the decay rate of the i-th particles including 
possible annihilation. 

According to the Smoluchowski approach, the reaction rate constant A;p Si R (let us 
take as an example reactions (1 1 6H1 7j) ) may be written as follows: 



/c PS)R (T,t) = 47r-D PSiR (T).Rp s>R 



, Rp s ,R 



(19) 



'irDp B) -B.t 

where i?p Sj R = -Rp s + R is the reaction radius of the i + Ps chemical reaction, i?j is 
approximately the geometric radius of z'-th reactant, Rp s ~ Ry + 1/x is the radius 
of derealization of the Ps wave function (1/x is the under-barrier penetration length 
of the Ps center-of-mass wave function; here, we consider Ps as a point particle in a 
potential well). 

The temperature dependence of the rate constant arises through the diffusion 
coefficients of the reagents, -Dp S]R (T) = Dp s +D^. According to the Einstein relationship 
and the Stokes formula, the diffusion coefficients may be expressed as: 

T T 

Dr = 6^(TJR> DPS = 4tc V (T)R Ps > (20) 

where r\ is the viscosity of the medium and R and Rp s are the hydrodynamic radii 
of R particle and Ps bubble. Although the applicability of these expressions on the 
atomic scale might be questionable, our studies in neat water [2U [23] have shown that 
it remains valid within a reasonable accuracy. Usually the viscosity rj(T) drops down 
approximately exponentially with T. 

The knowledge of the surface tension of the liquid together with the Ps bubble 
model allow one to estimate the equilibrium radius of the Ps bubble and calculate the 
pick-off annihilation rate X po (t,T), which is of utmost importance to interpret the PAL 
spectra: 

X po (t,T) = X p -P(Ru(t,T)). (21) 

Here, X p is the "free" positron annihilation rate and P(Ru) is the underbarrier 
penetration probability of e + into a bulk of the liquid (into space containing molecular 
electrons). Its calculation is discussed below. 

3.2. Low temperature region 

In this region, the Ps bubble formation time is larger than the e + and Ps annihilation 
lifetimes. On a scale of 1 ns, the size of a preexisting void in which a quasi-free Ps has 
been localized does not change, since Uubbie ^> 1 ns. This situation is similar to that 
existing in polymers. Typically at these T the diffusion time t<nff is much larger than 
o-Ps lifetime (~ 1 ns), so one may completely neglect the diffusion motion of radiolytic 
products. For example in glycerol it happens at T < 30 °C, Fig. [HJ 

To describe correctly the PAL spectra, one should average the pick-off annihilation 
kinetics over the size distribution of the preexisting voids in which Ps localization takes 
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place. This can be done by using the theory of free volume entropy fluctuations [2T] , 
which claims that every molecule possesses some free volume V{ with the probability 

ex exp(-Vi/v F ), v F (T) « v ws (T) - v vdW , (22) 

where Vp is a difference between the volume of the Wigner-Seitz cell and van-der-Waals 
volume of a molecule. Because of the approximate character of this expression, Vp may 
be considered as an adjustable parameter derivable by fitting the spectra in the low 
temperature region. 

If v is the volume of a preexisting void in which the Ps atom can localize, the 
probability to find such a void is proportional to exp(— v/vp). Neglecting the influence 
of any intratrack reactions (because of low-T and not discussing their possible tunneling 
nature), the calculation of the Ps pick-off annihilation kinetics reduces to averaging 
the exponents that relate to the pick-off annihilation of individual Ps atoms located at 
different voids: 

/ e -\ vo (Ru)t\ _ f°° e -X P o(Ru)t . e -v/v F m f& v= &KRu , 23 x 

Jv • v f' 3 

w u min *■ 

Here, v m i n is the volume of the minimal cavity which may trap Ps (section HJ). 

3.3. Intermediate temperatures 

The most complex task in the PAL data processing refers to the intermediate 
temperature region, where the bubble formation time is comparable with the e + lifetime 
in the medium. qf-Ps localizes in one of the preexisting voids, which further slightly 
increases in size; the time at which the equilibrium volume is reached depends on T. 
Obviously, the quantum-mechanical pressure exerted by Ps on the wall and thus, the 
growth of the void, cease when e + annihilates. 

At these temperatures, the Ps fate may also depend on the intratrack reactions. 
Therefore, the corresponding chemical kinetic equations must be solved by assuming 
that Ps initially localizes in a preexisting void of a given radius, which will increase with 
time until annihilation; the process must be averaged over the size distribution of the 
preexisting voids [8]. 

3.4- Interpretation the PAL spectra in liquid water at different T 

To solve the problem with a puzzling temperature behavior of the long-lived component 
in water, Fig. we proposed to reject the conventional exponential deconvolution of 
the spectra and explicitly take into account the intratrack processes mentioned above 
pH [28] . The decrease of r 3 vs. T was ascribed to the increasing efficiency of the oxidation 
( fT6l) and ortho-para-conversion ( fl7|) reactions between Ps and intratrack radicals (mainly 
OH-radicals). 

All parameters included in the model have a clear physical meaning: the Ps 
oxidation reaction radius R ox « Rp s + -Roh which enters the rate constant of the 
Ps oxidation reaction ffT6|) . relative contact density r\ c in the Ps bubble state, free 
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Figure 9. Left: Temperature dependence of the Ps oxidation reaction rate constant 
(at an infinite time, t — > oo) multiplied by the "initial" concentration of oxidizers 
c ox (0) in the e + blob (OH and H30 + , see HUES) f° r details). D ox is the sum of Dp s 
and the diffusion coefficient of the oxidizer (OH radical, -Doh ~ 1.5 • 10~ 5 cm 2 /s). 
The obtained dependence agrees well with the Stokes-Einstein law, ex T/r)(T). Right: 
T-dependences of the oxidation reaction radius, R ox ,Ps = Rox + ^Ps, the Ps radius, 
i?p s = Rjj + \j k and the relative contact density parameter r\ c in Ps. R ox is the radius 
of the oxidizer (approximately, the geometric radius of the OH radical, 1.4 A). 



positron annihilation rate, the Ps formation rate constant and some others (5 adjustable 
parameters in total). 

It is seen that the blob model explains well the experimental data in a wide range 
of temperatures and magnetic fields. It does not lead to contradictions with known 
radiation chemistry data. The agreement between theory and experiment became even 
better when taking into account the time dependence of the Ps reaction rate coefficients. 
The temperature dependence of the Ps diffusion controlled reaction rate constants agrees 
with the Stokes-Einstein law. As we have obtained. Good fitting of the PAL data is 
obtained when the Ps hydrodynamic radius is equal to the radius Rp s = Ru + 1/x of 
the Ps bubble (here 1/x rj 0.5 A is the underbarrier penetration depth of the Ps wave 
function) . 

4. Ps bubble models 



Typical lifetimes of a para-positronium atom in condensed medium are about 130-180 ps. 
They are close to the p-Ps lifetime in vacuum (125 ps). The ortho-positronium lifetime 
in a medium is considerably shorter (about 100 times; some ns) in comparison with that 
in vacuum. This is due to the so-called pick-off process - prompt 27-annihilation of the 
e + , composing Ps atom, with one of the nearest e~ of surrounding molecules, whose spin 
is antiparallel to the e + spin. Just this property turns Ps into a nanoscale structural 
probe of matter. The theoretical task consists in calculating the pick-off annihilation 
rate A po , i.e. in relating X po with such properties of the medium like surface tension, 
viscosity, external pressure and size of the Ps trap. 

Originally, to explain the unexpectedly long lifetime of the ortho-Ps atom in liquid 
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Figure 10. Tao-Eldrup model of Ps atom in a 
liquid phase. It is assumed that Ps is confined in 
an infinite spherically symmetric potential well 
of radius i?^. \l/(r) is the center-of-mass wave 
function of Ps. R is the free volume radius of 
the Ps bubble. 5 is the penetration depth of the 
molecular electrons into the Ps bubble. 



helium R.Ferrel [29] suggested that the Ps atom forms a nanobubble around itself. This 
is caused by a strong exchange repulsion between the o-Ps electron and electrons of the 
surrounding He atoms. Ferrel approximated this repulsion by a spherically symmetric 
potential barrier of radius Roq. To estimate the equilibrium radius of the Ps bubble 
he minimized the sum of the Ps energy in a spherically symmetric potential well, i.e. 
7r 2 /i 2 /4mi?^ = -if(iraB/Roo) 2 i Ry=13.6 eV, and the surface energy, AnR^cr, where 
a is the macroscopic surface tension coefficient. The following relationship is hereby 
obtained for the equilibrium radius of the bubble: 



TT 2 a 2 

g Ry + AttR^ct -b- min over R 
R 



Rr 



a B 



8aa 2 B 



i/i 



(24) 



4-1. The Tao-Eldrup model 

Ferret's idea got further development in the studies of Tao [30] and Eldrup et al. [31]. 
They considered the Ps atom as a point particle in a liquid, i.e. in a structureless 
continuum, Fig. (TDJ The repulsive Ps-liquid interaction was approximated by a 
rectangular infinitely deep spherically symmetric potential well of radius -R^. In such a 
well, the wave function of a point particle has the following standard expression: 

sin(7rr/i? 0O ) 



*(0 < r < R c 



mr > Rc 



0. 



(25) 



y/2-irRoo r 

Here, r is the Ps center-of-mass coordinate. Because the Ps wave function equals to zero 
at the bubble radius (and outside), there is no e + overlapping with outer electrons of a 
medium. So, pick-off annihilation is absent. To overcome this difficulty it was postulated 
that molecular electrons, which form a "wall" of the Ps bubble, may penetrate inside the 
potential well. This results in the appearance of a surface layer of thickness S = Roo — R 
having the same average electron density as in the bulk. As a result, the pick-off 
annihilation rate X po becomes non-zero. It is proportional to the e + overlapping integral 
with the electrons inside the bubble: 

/ \^(r)\ 2 4:7ir 2 dr = — 

Jr 



A 



po 



A+P 



R, 



P, 



Rr. 



(26) 



2vr 

This is the well-known Tao-Eldrup formula. Here, A + ~ 2 ns^ 1 is the e + annihilation rate 
in an unperturbed medium (it is proportional to Dirac's 27-annihilation cross-section 
and the number density of valence electrons). The thickness 5 of the electron layer is 
an empirical parameter, which may have different values in various media. 
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Figure 11. Dependence of the 
experimental pick-off annihila- 
tion rates [32, vs. surface tension 
in different liquids. The solid 
curve shows the correlation given 
by the Tao-Eldrup at A + = 2 
ns -1 and optimal value S = 1.66 
A (obtained from fitting of these 
data by means of Eq. (|2"6"|) ). The 
dashed curve illustrates the sim- 
plest approximation X po ex a 1 ' 2 . 



Substituting Eq. (El) for R^ into Eq. (T26 
X po and o with one adjustable parameter, 5. 



, one obtains the relationship between 
It may be easily obtained by fitting 
experimental pick-off annihilation rates with the relationship (126]) . Fig. [TTJ Thus we 
obtain § rj 1.66 A. Eq. ( |26l) with this value of S is widely used for recalculation of the 
observed pick-off annihilation rate into the free volume 4ttR 3 /3 of the cavity, where Ps 
atom resides and annihilates. 



4-2. Further development of the Ps bubble models. "Non-point" positronium 

Along with the development of the "infinite potential well" Ps bubble model, another 
approach based on the finite potential well approximation was also elaborated [33j EH 
|35| |36| [32] . However in both approaches, the Ps atom was approximated by a point 
particle. This leads to a significant simplification, but it is not justified from a physical 
viewpoint, because: 

1) the size of the localized state of Ps (size of the Ps bubble) does not significantly 
exceed the distance between e + and e~ in Ps; 

2) during the formation of the Ps bubble there is a substantial variation of the 
Ps internal energy (particularly of the Coulombic attraction of e + and e~), which is 
completely ignored in the "point-like" Ps models. In a vacuum or in a large bubble, the 
internal energy of Ps tends to — Ry/2 = —6.8 eV. In a continuous liquid (no bubble) 
with the high-frequency dielectric permittivity e ~ n 2 (n rj 2-3 is the refractive index) 
the energy of the Coulombic attraction between e + and e~ decreases in absolute value 
by a factor e 2 ~ 4-9. The same takes place with the total Ps binding energy, which 
tends to the value — Ry/2e 2 ps —(1-1.7) eV (this is a simple consequence of the scaling 
e 2 — > e 2 /e of the Schrodinger equation for Ps atom). Thus, the change in the Ps internal 
energy during Ps formation may reach 5 eV. Obviously, this represents an important 
contribution to the energetics of Ps formation. The aim of the present work is towards 
a more accurate estimation of this contribution, that has not been done yet. 
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There is only a small number of papers where the consequences of the finite size 
of Ps are discussed in application to positron annihilation spectroscopy. To calculate 
\ po , the Kolkata group [37] suggested to smear the Ps atom over the relative e + -e~ 
coordinate exactly in the same way as it is in a vacuum. Such an approach is valid 
for rather large bubbles. However, they do not discuss the variation of the internal Ps 
energy. 

In [38] the Ps atom is considered as a finite sized e~e + pair, but the variation of the 
Coulombic interaction because of dielectric screening is not discussed. It was assumed 
that e~ is confined in an infinite potential well and e + is bound to it by means of the 
Coulombic attraction. The wave function of the pair was taken as a series of orthogonal 
polynomials, their weights being determined from a minimization procedure of the total 
energy of the pair. 

4-3. Hamiltonian of e + e~ pair in a medium 

Let the e + e~ pair (Ps atom) have already formed in a liquid a nanobubble (spherical 
cavity; Ps bubble) of radius R (the onset of coordinates is taken at the center of the 
bubble, Fig. I4.3p . Together with the molecules surrounding the e + e~ pair, one has to 
deal with a quite intricate many-body problem with a complex hamiltonian. We reduce 
it to the following form: 

H « - /l2(A 2 + J A ' ) + U(r+) + U{vJ) - U c (r + , r , R, e). (27) 

Terms with Laplacians A + and A_ over r + and r_ (e + and e~ 
coordinates) stand for the kinetic energies of the particles. U(r + ) 
and U(r_) describe the individual interaction of e + and e~ with 
the medium. For them we adopt the following approximation: 

TTt \ / ' r + < ^' TTt \ I ° ' T - < R > foa\ 

U{T+) = J K, r + > R, "^ = { V4-, r_ > R. < 28) 

Here, V^~ and V$ are the e + and e~ work functions, respectively (Vo is a commoner 
notation for the electron work function). The work function is usually introduced as 
the energy needed for an excess particle to enter the liquid without any rearrangement 
of its molecules and to stay there in a delocalized state, having no preferential location 
in a bulk. One may say that Vq and V _ are the ground state energies of the quasifree 
e + and e~, because their energies at rest after having been removed from the liquid to 
infinity are defined to be zero. 

Vq~ consists of 1) the e~ kinetic energy, arising from its exchange repulsion from 
the "core" electrons of molecules (atoms), and 2) the energy due to the polarization 
interaction of e~ with the mediumjjj According to the theory of the quasifree electron 

X In case of e + the kinetic contribution to V + is due to the Coulombic repulsion from the nuclei (the 
exchange repulsion is absent). 
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Table 2. Electron work function for different liquids at room temperature [5] 



Liquid 


V -, eV 


Liquid 


V ~, eV 


helium; 4.2 K 


1.3 


benzene 


-0.14 


n-dodecane 


0.2 


isooctane 


-0.17 


n-decane 


0.18 


toluene 


-0.22 


n-heptane 


0.12 


neopentane 


-0.38 


n-hexane 


0.1 


MeOH, EtOH, PrOH 


-0.4 


nitrogen; 77.3 K 


0.05 


xenon; 170 K 


-0.57 


n-pentane, c-hexane 


0.01 


water 


-1.2 


argon; 86.4 K 










approximately the same as for e 
eV. Thus |V + + V - 



[39], this polarization interaction may be estimated as a sum of two parts: a) interaction 
of the e~ with the molecule where it resides, U 1 ™ 1 (to calculate U™* the electron is 
considered as an electron cloud smeared over the molecule), and b) interaction of the 
e~ with all the other molecules, t/f"* = (1 — l/e)e 2 /2Rws, (this expression is similar to 
the well-known Born formula for the electron solvation energy). 

Experimental values for V ~~ are known for many liquids (Table Wj- Because of 
a lack of experimental data on the e + work functions, we shall admit that they are 

g + pa V Q . So we may conclude that \V + + V ~\ < 1 
is less than the variation of the internal energy of the pair, 
~ Ry(l — l/e 2 )/2 « 5 eV, related with the variation in the dielectric screening of 
the e + -e~ attraction in the bubble formation process. 

Note that the use of Eqs. ([28"]) for the potential energies of the e + and e~ interaction 
with the medium, assumes that the polarization interaction remains the same whether 
e + and e~ (both in the quasi-free states) are well separated or form the quasi-free 
Ps atom. Since, for distances larger than the size of a molecule, qf-Ps is nearly an 
electrically neutral particle, the contributions U™ 1 ~ U™ 1 , which come from a long- 
range polarization interaction of the quasi-free e + and e~ with the medium, should 
be absent in U(rJ) + U(r + ) in Eq. ( J27l) . Therefore, it is reasonable to consider at 
least two cases: 1) when the above mentioned polarization correction is neglected and 



Vr, + V+ — > and 2) when the terms U. 



OU I 



U 



out 



— 1 eV are subtracted from the 



work functions and V n 



Vq — > 2 eV. Both cases are considered below. 



In Eq. (127)1 U c stands for the Coulombic interaction between e + and e~ in a 
polarizable medium. Assuming that the medium has the dielectric permittivity e of 
the bulk and a spherical cavity of radius R (inside the cavity e — 1), one may calculate 
U c by solving the Poisson equation. Denoting the e + and e~ coordinates as r + and 
r_, U c may be written in the form of the following series via the Legendre polynomials 
P / (x = cos^) (HI: 



U c (r + < R, r_ < R) 2a B 
Ry r 



- 1 



2ob 
~R 



1 



oc 



l + l)Pi{x) 
^ 1 + l + l/e 
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Figure 12. Dependence of the e + -e~ Coulombic interaction energy for different 
locations of e + and e~ around the bubble (here we adopt the radius of the bubble 
R equal to 4as ~ 2 A). z + and z_ are the e + and e~ displacements from the center 
of the bubble along the z-axis. The dashed curves describe the unscreened (red) and 
completely screened (e times less; blue) Coulombic energies between e + and e - . 
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Here, the argument of the Legendre polynomials is x = cosf?, where t? is the angle 
between the z axis and the direction of r_. Note that the summation of these series is 
simplified considerably when using the following recurrent relationship 

P t (x) = [(21 - l)zft_i(z) - ( l ~ l)Pi-2(x)}/l. 

Particular dependencies of U c for some selected arrangements of e + and e~ and the 
cavity are shown in Fig. CGJ Thus, we are able to take into account the variation of the 
e + e~ Coulombic energy during the formation of the Ps bubble. Similarly, the dielectric 
screening is used in the polaron problem and the ion-electron recombination problem 
(Onsager's formula) [¥Ij W2\. 

4-4- Wave function of the e + e~ pair and minimization of its total energy (H) 

Keeping in mind further use of the variational procedure, let us choose the normalized 
e + e~ wave function in the following simplest form: 

exp(-r/2a - r cm /2b) r+ + r_ 



*+-(r+,r_; 



-r_.(29) 



In both cases of a rather large bubble and a uniform dielectric continuum, \l/+_ breaks 
into a product of two terms: the first one depends on the distance r between e + and e~, 
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and the second one depends on the center-of-mass coordinate r cm . Parameters a and b 
are the variational ones, over which we have minimized the energy of the e + e~ pair: 

E(a,b,R) = (^+_|#|^+_) ^min => a(R), b(R). (30) 

The simplest verification of the calculations is to recover two limiting cases. In case 
of large bubbles (R — > oo), one should reproduce the "vacuum" state of the Ps atom: 
its total energy must tend to — Ry/2 = —6.8 eV, the kinetic energy to +Ry/2 and the 
Coulombic energy to — Ry. In case of small bubbles (R — > 0), the delocalized qf-Ps 
state must be reproduced. The Schrodinger equation for qf-Ps has the same form as 
for the vacuum Ps, but with the substitution e 2 — > e 2 /e. Then the total qf-Ps energy 
tends to V^ + V ~ - Ry/ 2 ^ 2 , its kinetic part tends to +Ry/2e 2 = 1.7 eV (e = 2) and 
the Coulombic energy tends to — Ry/e 2 = —3.4 eV. Fig. [13] displays optimal values 
of a and b as well as different contributions to the total energy of the e + e~ pair when 
V+ + Vq = and 2 eV. 

4-5. Relative contact density and pick- off annihilation rate 

Using the wave function ( J29l) it is easy to obtain the relative contact density n c in the 
Ps atom: 



j j d 3 r + d 3 r_\^ + _(r + ,rJ)\ 2 5(r + — r_] 



a' 



:-! 



/c J/d 3 r + rf 3 r_|^^(r+,r_)| 2 5(r + -r_) a 3 (R)' K ' 

This quantity determines the observable Ps annihilation rate constant (including the 
case with applied permanent magnetic field). The resulting dependencies of rj c are 
shown in Fig. [13] (on the left). Because, for qf-Ps, parameter a is equal to £a#, for qf-Ps 
the value of r] c should be 1/e 3 = 1/8, which is well recovered in numerical calculations. 
When R increases, r] c approaches unity, because a tends to its vacuum value as- 

Knowing the expression for the wave function ff29l . one may calculate the positron 
overlapping Pr with molecular electrons, surrounding the Ps atom, and therefore find 
out the pick-off annihilation rate constant: 

X po (R) « X + P R , P R « / d 3 r + [d 3 r_ |tf+_(r+, r_)| 2 . (32) 

Jr+>R J 

Here, A + ~ 2 ns _1 is the annihilation rate constant of "free" positrons. Results of 
calculations of X po (R) for optimal a and b values, which correspond to the minimal Ps 
energy at a given R, are shown in Fig. [141 

4-6. Discussion on the non-point-Ps approximation 

It is usually considered that Ps is a solvophobic particle, i.e., it forms a bubble when 
entering a liquid because of exchange repulsion between e~ in Ps and the surrounding 
molecular electrons. If the work functions of e + and e~ are negative (Vq~ ~ Vq < 0), 
each particle considers a cavity as a potential barrier. So they are pulled to the bulk 
by polarization interaction with the medium. Nevertheless, even in this case the Ps 
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Figure 13. Dependencies of the optimal parameters a and b vs. R, the bubble radius. 
They enter the e + e _ wave function and yield the minimum of the total energy (H). 
The relative contact density r\ c and different energy contributions to (H) (at optimal 



a and b) are shown as well. The upper drawings correspond to the case Vq + V 







and the lower ones to Vn 



v n - 



2 eV. In both cases it was assumed that e 



bubble may be formed due to an enhancement of the Coulombic e + e~ attraction inside 
the cavity (no dielectric screening inside). This feature cannot be taken into account 
when Ps is simulated as a point particle. 

It is seen that the behavior of the total energy of the pair (red curves in Fig. [T3|) 
strongly differs from the Tao-Eldrup prediction (green dashed curves; the first term in 
Eq. (|2^|) . where i?^ is replaced by R), as well as from the expectation based on the 
finite potential well model (brown curves in Fig. [131 the Coulombic potential cannot be 
approximated well by a rectangular spherically symmetric potential). The same is true 
for the pick-off annihilation rate, Fig. [TU 

Calculations demonstrate one common feature: up to R < 1.5 — 2.2 A all 
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Figure 14. Pick-off annihilation 
rate constant of Ps, localized in a 
bubble of R when V+ + V- =0 
(green curve) and V + + V- = 2 eV 
(brown curve). For small R (< 2 
A) the calculated values of \ po are 
equal to A+. The red line shows 
pick-off annihilation rate constant, 
calculated according to the Tao- 
Eldrup formula. The dashed line 
is the calculation according to the 
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comparison we adopted that the 
depth of the well is (l-l/e 2 )Ry/2 « 
5.1 eV and its radius is R. The 
minimal radius of the well when 
there appears an energy level is 
Su =0.96 A). 



dependencies remain the same as in a medium without any cavity, but at larger R 
there are significant deviations. This is related to the known quantum mechanical 
phenomenon - absence of a bound state of a particle in a small finite 3d-potential well. 
In such cavities, Ps cannot be bound, it does not exert any repulsive pressure on their 
walls and does not stimulate their transformation towards the equilibrium Ps bubble. 
The possibility of finding a suitable preexisting cavity, sufficient at least for preliminary 
localization of qf-Ps, may be a limiting factor for the formation of the Ps bubble state. 
One may find an equilibrium Ps bubble radius by minimizing the sum of the total 
e + e~ energy (H) and the surface energy of the bubble. For water it turns out to be 5-5.2 
A which is about 2 A larger than predicted by the Tao-Eldrup model. For such a large 
bubble, the relative contact density is r] c « 0.9, Fig. [131 It is somewhat higher than 
the experimental values (0.65-0.75 [28]). This discrepancy may indicate that e + and e~ 
really interact with a medium in a different way, for example, Vq > Vq~. It means that 
the Ps electron may be trapped by a cavity, and e + will be bound to this trapped e~ 
by the Coulombic attraction. This scenario may be also considered in the framework 
of the present approach, but the expression for the trial wave function of the new pair 
must be written in an "asymmetric" (towards e + and e~) form: 

]/2a-r^/2b) 



*_ 



exp( 



r+, r_) « ^" v r+ ^T^r •-'-' . (33) 

87rva 3 & 3 

Any Ps bubble model reduces the original many-body (multi-particle) problem 

to a simpler one, that of one or two particles in an external field, which simulates 

the interaction with the medium. To calculate this field one usually relies on some 

macroscopic approaches. However, their validity always remains uncertain (for example, 

how to relate the actual arrangement of molecules around the Ps bubble with the jump 

of dielectric permittivity outside the bubble and so on). 
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5. Conclusion 

In this chapter an attempt has been made to trace the fate of a positron in molecular 
liquids, starting from its entering the medium till annihilation with an electron. 

Energetic restrictions have been shown to leave only one possibility to form Ps, 
when the Ps atom is formed as a result of combination of the thermalized positron with 
one of the track electrons (some it is called as the (re) combination mechanism or the 
intratrack mechanism). This reaction proceeds in the e + blob, terminal part of the 
e + track. Considering the ionization slowing down of the projectile positron, one may 
conclude that the e + blob contains several tens of ion-electron pairs, which are confined 
in a spherical region with a radius of several nanometers. Ion-electron pairs in the blob 
are bound together because of their mutual Coulombic attraction. This is why further 
expansion (out-diffusion) of the blob is characterized as the ambipolar diffusion. 

When the positron energy becomes less than ^Ry=13.6 eV (subionizing e + ), this 
particle can escape from the blob. It becomes thermalized either outside the blob or 
inside it. Being inside the blob it may pick up intratrack electron and form Ps. 

Primary ion-electron pairs transform into secondary radiolytic products (typical 
concentration > 0.01 M), which are chemically active and thus can interact with Ps: 
oxidize it and/or stimulate its ortho-para conversion. Because the number of ion-electron 
pairs is large and the diffusion displacement of the species on a nanosecond timescale is 
comparable with the size of the blob (at least at a high-T region), it appears inescapable 
to describe these processes on the basis of inhomogeneous diffusion kinetic equations in 
terms of concentrations of the species, rate constants, diffusion coefficients etc. 

In the case of water it is clearly seen that taking into account interactions with 
intratrack products (Ps oxidation by OH radicals and H30 + ) is important. Otherwise, 
it is not possible to explain the experimentally observed decrease of the ortho-Ps lifetime 
with temperature. 

Calculating the pick-off annihilation rate constant is of the utmost importance. 
This is the main goal of the Ps bubble models. Up to now most of calculations deal 
with the "point-like" -positronium, by using its center-of-mass wave function only. In the 
framework of this approximation it is not possible to take into account the variation of 
the internal (Coulombic) energy of the e + -e~ pair during Ps bubble formation. However, 
we have shown that this energy contribution is not small and may play a decisive role. 

Finally, we would like to draw attention to one of the possible applications of the 
positron spectroscopy: fast detection of potentially carcinogenic chemical compounds. 
These substances are considered as one of the major causes of human cancer. Modern 
technologies have led to new potential chemical carcinogens with which people may 
be in contact in everyday life. The European Community considers that about 100 
such compounds need to be screened every day. Therefore, there is a need for a fast 
and cheap method of testing such compounds. Positron annihilation spectroscopy may 
provide such a method [431144]. 
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